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Abstract 

I give a proof, using the unitarity-based method, of the collinear factorization of the 
leading-color contribution to gauge-theory amplitudes. The proof also provides a concrete 
formula which can be used to compute the associated splitting amplitudes. 
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1. Introduction 



Gauge-theory amplitudes have simple and universal features in various singular limits. In the 
collinear limit, when the momenta of two massless particles become proportional (k a — > z(k a + k b ), 
kb — > (1 — z)(k a + fc&)), we may summarize these factorization properties in a set of universal 
splitting amplitudes. The squares of the tree- level collinear splitting amplitudes, summed over 
helicities, yield the Altarelli-Parisi kernels [1,2]. Their integrals play an important role in the 
infrared cancellations that are needed for calculations of physical quantities at next-todeading 
order in perturbative quantum chromo dynamics [3,4,5,6]. 

Beyond leading order, splitting amplitudes also have relatively simple forms. Their explicit 
forms at one loop [7,8] were extracted by taking the collinear limit in various five-point one-loop 
amplitudes. They are useful as complements to the unitarity-based method of computing one-loop 
amplitudes [9,10], and will play a role in the cancellation of infrared divergences in next-to-next- 
to-leading order calculations. 

Bern and Chalmers [11] have previously given a proof of collinear factorization at one loop. 
It has the virtue of being quite general, and not limited to gauge theories. It is unclear, however, 
how to generalize their proof beyond one loop. Furthermore, even at one loop it provides neither a 
simple understanding of the structure of the splitting amplitudes, nor a simple method of computing 
them. The purpose of this paper is to provide an independent proof of collinear factorization which 
generalizes readily to higher loops. It will also provide a simple and concrete formula for computing 
the splitting amplitudes at one and two loops. The proof is an application of the unitarity-based 
method for computing loop amplitudes developed by Bern, Dixon, Dunbar, and the author [9,10]. 

2. Order of Limits 

The properties of non-Abelian gauge-theory amplitudes in the collinear limit are easiest to 
express in the context of a color decomposition [12]. For tree- level all-gluon amplitudes in an 
SU(N) gauge theory the color decomposition has the form, 

A t : ee ({k i ,X i ,a i })= Tv(T a ^ ■■■T a ^)A t ^ ee (a(l x \...,n x ")), (2.1) 

ae.s n /z n 

where S n /Z n is the group of non-cyclic permutations on n symbols, and j x i denotes the j-th 
momentum and helicity. As is by now standard, I use the normalization Ti{T a T b ) = 5 ab , and I 
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suppress powers of the coupling constant. One can write analogous formulae for amplitudes with 
quark-antiquark pairs or uncolored external lines. The color-ordered or partial amplitude A n is 
gauge invariant. 




Figure 1. The discontinuity of a loop amplitude expressed as a product of tree amplitudes. 
In the cut-based approach to one-loop calculations, we have the following basic equation, 

,.1-loop _ f d 4 2£ l i tree [ .tree 

2^ / ( 2 7r) 4 - 2e £ 2 left U-KY^ ( 2 - 2 > 

cuts K 2 J y ' y ' 

where the loop integral is promoted from the integral over phase space, and the sum over all cuts 
is understood to count only once integral functions that have cuts in multiple channels. The right- 
hand side is depicted in fig. 1. (I will leave the sum over all intermediate particle types and physical 
polarizations implicit in the cut expression. In the four-dimensional helicity scheme [13,14], these 
would be the positive and negative helicity states; in the conventional dimensional regularization 
scheme [15], it would include 'e' helicities as well.) 

In general, we are interested only in the resulting terms through finite order in the dimensional 
regularization parameter e. Reconstructing the amplitude by summing over the cuts does not 
however commute with expanding around e = 0; only if all cuts satsify a power-counting criterion 
can the full amplitude, to order e°, be reconstructed from the D = 4 tree amplitudes on either side 
of the cut. Otherwise, we must keep the cuts to higher orders in e (or use other tricks) in order to 
obtain the finite rational terms in the amplitude. 

The singular limit in which we are interested poses additional interchange-of-limits problems. 
The limit in which k a || fe&, or more generally in which a single invariant s a b = (k a + kb) 2 of two 
color-adjacent legs vanishes while all other invariants stay (roughly) fixed, fails to commute with 
the expansion around e = 0, as can be seen in the example, 



(As is conventional in investigating infrared divergences, I will take e < 0.) This interchange-of- 
limits issue in specific integral functions is addressed by the Bern-Chalmers discontinuity functions 
[11]. 

In what order should we take the singular and four-dimensional limits? It might seem more 
natural to take the e — > limit first, but this is not the best choice. The problem is precisely that in 
this order, we lose control not only of divergent terms but also of finite terms in the integrals (the 
In 2 s above). Keeping them under control would mean keeping track of all finite terms in integrals, 
a rather formidable task. On the other hand, if we take the singular limit first, then all one-loop 
integrals with more than one external invariant (i.e. all but one-mass triangles and bubbles) are in 
fact finite in any singular limit. This may seem bizarre — what about the st in the denominator 
of the one-mass box, for example? — but is easily seen to be true using the Feynman parameter 
representation. (The resulting collinear-limit integral may, of course, be more singular in the later 
e — ► limit, but that is a later concern.) Indeed, it is instructive to consider the s, t — > limit of 
the reduced one-mass box [16]; what naively looks to be uncontrolled growth the standard e — > 
expansion, 




(with rr = T(l + e)r 2 (l — e)/T(l — 2e)) turns out to be quite well behaved (after analytic contin- 
uation) if we go back to the Feynman parameter representation, 



d\ 



(^r(2 + e) f d^a - ^ E ^ } 2 We — (^ 2 )T(2 + e) f 

Jo [— sai<23 — ta2CL4 — TrrcuaiJ Jo [— m^a^ai] 

= -2rrV) £ (-m 2 )- 2 - £ 

(2.5) 

The analytic continuation involved, however, is not appropriate for us, as it amounts to subtracting 
poles in s and t. (The scale \i is introduced by dimensional regularization and would eventually 
become the MS renormalization scale.) 

We do want to pick up such poles, so the appropriate limit for our purposes is given by 
expanding the integrated form, leaving epsilonic powers alone, and taking care to drop only terms 
suppressed by integer powers of s/m 2 or t/m 2 . 



Using the form of l\ m expressed in terms of hypergeometric functions, 
2r r 



e 2 st 



st \ e / s + i - m 2 

fj, 2 (t — m 2 )J 2 \ ' ' ' t — m 2 

st \ e „ / „ s + t — to 2 



+ — ^ 2F1 -e,-e;l-e; 2 (2.6) 

sim 2 \ 6 „ / (s + i — m 2 )m z 

' 2 Fi -e,-e;l-e; 



(s + t — m 2 )m 2 — st J \ ' ' ' /U 2 [(s + t — m 2 )m 2 — st] 

dropping terms suppressed by integer powers of (say) s/m 2 , and only then expanding in e, we find 

£ (^) ^ 1 " " " " ¥ {^h) " r(1 " f)r(1 + f) 

e 2 \ —m z ij, z J 3 

which is well-behaved in the singular limit. We should thus take the singular limit first, and only 
then the four-dimensional one. In the following discussion, unless otherwise stated, I will then keep 
amplitudes and cut expressions to all orders in e, unexpanded in a series. 

Of course, there is no free lunch. The discrepancy between eqn. (2.5) and (2.7) hints at the 
problem that naive manipulations with the Feynman-parameter representation dispose of terms 
like s~ e , which we really want to keep. If we are considering a cut in one channel, say s±, and 
the resulting cut is proportional to s^ e , where S2 — > in the desired singular limit, then it would 
appear at first glance that this contribution will disappear quietly, and we won't pick it up. (Of 
course, if the s^~ e does not show up as part of an integral function containing a cut in the S\ 
channel, but merely as a separate function from the integral reductions, then it should be dropped 
anyway because the cut in the si channel necessarily has an ambiguity of that form.) However, 
this contribution will also have a cut in the S2 channel, and we can expect to pick it up there. We 
will therefore have to be very careful about taking the limit of the cut in the singular channel. This 
subtlety would also mean in considering singular limits with more than one singular invariant, that 
we would have to take cuts in multiple channels simultaneously, in order to pick up contributions 
proportional to sps^ e . 

While we cannot take the singular limit before computing the cut for the cut in the singular 
channel itself, what about other channels? In general, an inverse integer power of the singular 
invariant must be supplied by the amplitudes on either side of the cut, or else by the loop integral 
whose discontinuity is given by the cut. (The helicity algebra eventually softens this singularity 
to a square-root from a full pole.) If the inverse integer power is supplied by an amplitude on 
either side of the cut, we can clearly take the singular limit before considering the cut. The same 



turns out to be true (except in the singular channel, as discussed above) even if the inverse integer 
power arises from the loop integral itself. Because we are working to all orders in e, the integral 
will necessarily have a cut in the singular channel. Thus whatever mistake we would be making 
by taking the limit naively in some other channel will be fixed up when we consider the cut in 
the singular channel itself. This implies, in particular, that cuts in which the collinear momenta 
are on opposite sides of the cut can be dropped. Such configurations do give contributions in the 
collinear limit, for example in the form of the hard two-mass box integral; but this contribution 
can be picked up in the singular channel. 



3. Collinear Factorization at One Loop 

At one loop, all amplitudes can be written as sums of leading-color amplitudes*. It is thus 
sufficient to consider the collinear limit of leading-color amplitudes in order to characterize the 
complete collinear behavior of a one-loop amplitude. Furthermore, the amplitude is finite in the 
limit unless the two collinear legs are (color-) adjacent, so restrict attention to this case. 

The amplitude can be reconstructed completely from the knowledge of all cuts; because we are 
working to all orders in e, there is no possible rational ambiguity. If we can determine the collinear 
limits of all cuts, we will then have determined the collinear limit of the amplitude. 

Consider first a generic cut, in which the two collinear legs a and b are a proper subset of the 
legs on one side of the cut, and label the legs on the other side from c through d. The cut in this 
channel, t c ,,,a = (k c H + kd) 2 , is then 

A^(l,...,a,b,...,n)\ tc d cut = 

J cOPS 4 " 2 ^!, -£ 2 ) A^ m+2 (£i, c, ...,d, -£ 2 )Al e ; 2 (e 2 , d + l,...,a,b,...,c-l, -£,) 

(3.1) 

To obtain the collinear limits of this expression, we must make use of the tree-level factorization, 
< ee (l, . . . , a A » , &\ . . . , n) k -^° Yl SpUt^(a A - , 6^)^(1, . . . , (a + by, . . . , n) , (3 2) 

ph. pol. a 

where Split tree is the usual tree splitting amplitude, the notation 'a + V means k a + kb, and where 

, 1/2 

the notation 'ph. pol.' indicates a sum over physical polarizations only. Split ree (a, b) goes as s ab ; 
terms not singular in the limit are omitted in this factorization. This factorization is depicted 
schematically in fig. 1. 

* Or primitive amplitudes, when some of the external legs are fermions. 
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This form may be derived from either the Berends-Giele recurrence relations [17], or else [18] 
from the Koba-Nielsen open-string amplitude [19]. Either derivation shows that this factorization 
holds only for on-shell (that is, physically polarized) legs a, b, but in arbitrary dimension. The 
following arguments will thus go through equally well in the four-dimensional helicity scheme, the 
conventional dimensional regularization scheme, or the original 't Hooft-Veltman scheme. 




a\\b 



a 




a+b 




Figure 1. A schematic depiction of the collinear factorization of tree- level amplitudes, with 
the amplitudes labelled clockwise. 



The tree splitting amplitude is given by the appropriate limit of the three-point Berends-Giele 
current, 



Splht ee (a A * , b x ») = -J— 4 A «) • e { b Xb) (k b - k a ) ■ e W + 2k a • 4 Ab) 4 A «) • - 2k b • e { b Xb) ■ e 



V2~S ab 

where S denotes the fused leg, ks = k a + 
In the limit, eqn. (3.1) then yields, 



(3.3) 



E 

ph. pol. a 



Split^ 6 (a A » , b Xb ) J dLIPS 4 " 2e (h , 



x Cm+2(<i> c, ...,d, -e^A^s^, d+i, . . . , (a + by, . . . , C-l, -h) 

2 Split*? (a x ",b x ») A 1 ^(l,...,(a + by,...,n) 

ph. pol. a 



(3.4) 



t c ... d cut 



As noted in section 2, we need not consider cuts where the momenta are on opposite sides of 
the cut (in which case they are both necessarily adjacent to it). The above derivation breaks down, 
as expected, if a and b are the only legs on one side of the cut; but all contributions except those 
detectable in the singular channel take the form presented in eqn. (3.4). This leaves us with the 
singular channel, which I consider next. 
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4. The Singular Channel 



The cut in the singular channel has the following form, 
A^(l,...,a,b,...,n)\ sab cut = 

£ J oTiPs 4 - 2e (^ 2 )< e ^r^+i,...,a-^ 

ph. pol.Ai 2 

(4.1) 

(I have relabelled l\ —>■ —i\ in this equation.) We already know we must treat it carefully, in 
particular we are not allowed to take the singular limit before performing the cut integral. It is 
nonetheless reassuring to see that the expression itself raises several warning flags, independent of 
the discussion in section 2. The cut integral is somewhat peculiar, because the volume of phase 
space vanishes in the singular limit. Furthermore, the factorization we used in the previous section 
fails. We cannot factorize the four-point amplitude into splitting amplitudes times three-point 
amplitudes, because the three-point amplitude vanishes on-shell. (Note, incidentally, that this does 
not mean that the four-point amplitude is not singular in this limit; on the contrary, it is more 
singular, having a full pole rather than just a square-root singularity. This is of course just the 
usual forward-scattering singularity of a gauge theory.) 

The other special feature of this cut is what allows us to simplify this expression. In this cut, 
momentum conservation conservation forces S£ 1 e 2 — > in the limit. This in turn allows us to use 
the tree- level factorization (3.2) on the amplitude on the left-hand side of the cut (4.1), 

A^ e (-l-^ , b + 1, . . . , a - 1, -£~ x *) 

£ spii^r (-e^, -r^) A^a-ii - i 2 y,b + 1, . . . , a - 1) 

ph. pol. a \ * / 

= Split*™ Mj Aa , ~^ Xl ) + bY, b + 1, . . . , a - 1) . 

ph. pol. a 

Note that the only remaining dependence on the cut momenta is in the tree-level splitting function. 




The cut in eqn. (4.1) thus becomes 



£ + a-1) 



ph. pol. 



J dLlPS 4 - 2e (£ 1 ,e 2 ) Split*^ (-^2 Aa ,-^r Al ) A r e (^2 a » aAo » &Ai »^i 1 ) ( 4 - 3 ) 



ph. pol. IT 



Sab CUt 



in the collinear limit. Since this channel is the only one in which Split ~ oop has a cut, we can 
immediately write down a formula for it, 



Split it° op (a Ao ,6 Ai ) = 

E / (gS h + « + b r\-r^ (e + ,; + kb? (4.4) 

ph. pol. Ai,A 2 J y ' y 1 

x Al ree ((-t-a-b) X2 ,a x «,b Xb ,e x i) 



The restriction to physical polarizations is important; it will give rise to transverse projection 
operators inside the loop. The defining equation is depicted graphically in fig. 3. 

If we now combine all cuts, we find in agreement with the known result [7], also shown in fig. 

2, 



^ loop (l,...,a A »,5 A6 



£ (Splits (« Aa , b Xb ) 4-T P (l, • • • , (a + by, . . . , n) 



(4.5) 

+ SplrC OOI V" , • • • , (a + &r, • • • , n) v 



ph. pol. cr 

,1-loop/ A a iA \ /itree 



This completes the proof of the universality of collinear factorization at one loop. Unlike the 
proof of Bern and Chalmers, it provides a compact formula for computing the one-loop splitting 
amplitudes. (Explicit computations with this formula will be presented elsewhere [20].) As we shall 
see in following sections, it also generalizes nicely beyond one-loop amplitudes. 
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Figure 2. A schematic depiction of the collinear factorization of one-loop amplitudes. 
5. Two-Loop Collinear Factorization 

The analysis in the previous sections extends readily to two- loop amplitudes. At two loops 
or beyond, leading-color amplitudes (or equivalently, planar color-ordered amplitudes) do not suf- 
fice to construct the full amplitude. In the following, I will consider explicitly only leading-color 
amplitudes, though the analysis can be extended to subleading-color ones. 

The basic equation (2.2) of the unitarity-based method has an analog at two- (or higher-) 
loop order. The sum over cuts now includes not only two-particle cuts but three-particle cuts as 
well. The amplitudes on either side of the cut are no longer just tree amplitudes, but may be loop 
amplitudes as well, 

Al-^(l,...,a,b,...,n)\ tc d cut = 

J dLIPS 4 - 2 ^,^)^-^^ 

+ J dLIPS 4 " 2 ^ , £ 2 ) <! e m+2 (^i, c, . . . , d, £ 2 )A 1 m l ^(-£ 2 ,d + l,...,a,b,...,c-l,-£ 1 ) 

+ j dUPS 4 - 2 ^e 1 ,£ 2 ,£s)A t :i e m+3 (h,c,...,d,£ 2 ,£s) 

x Al%(-£ 3 , -£ 2 ,d + 1, . . . , a, b, . . . , c - 1, -h) . 

(5.1) 

Note that in the three-particle cuts, all the £i must be adjacent in order to obtain a leading-color 
(planar) contribution; in the labelling used here, £3 is actually the middle leg crossing the cut. 
The sum over intermediate particle types and physical polarizations should again be understood 
implicitly. 

The arguments on the ordering of limits given in section 2 are in fact quite general, and apply 
to higher-loop amplitudes as well. We may thus take the singular limit first, in all channels other 
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t c ... d cut 



, n , 



t c ... d cut 



(5.2) 



than that of the singular invariant s a &. To do so, we must use both the tree-level collinear limit (3.2) 
as well as the one-loop limit (4.5), whereupon we find that the above expression has the following 
limit, 

Y Spht^V^) JdUPS 4 -^(e 1 ,£ 2 )A 1 r ; l _z p + 2(h,c,...,d,£ 2 ) 

ph. pol. a 

x^ e + e 1 (-^,d+l,...,(a + 6r,...,c-l,-^ 1 ) 

+ J dLiPs 4 - 2 ^£ 1 ,e 2 )A t :i e m+2 (e 1 , c ,...,d,e 2 ) 

xA 1 m l °° p (-e 2 ,d + i,...,(a + by,...,c-i,-e 1 ) 
+ j dups 4 - 2 ^e 1 ,t 2 ,£ 3 )A t :i e m+3 (e 1 ,c,...,d,e 2 ,t 3 ) 

x 4T+ 2 (-4, -*2,d + i,...,(a + by,...,c-i,-e 1 ) 

+ Y, Spht^° P (a A %^) |^LIPS 4 - 2e (^,^ 2 )<T m+2 (4,c,...,d,£ 2 ) 

ph. pol. a 

xA t ^ 1 (-e 2 ,d + l,...,(a + by,...,c-l,-£ 1 ) 
= Y Split^ 6 ^,^) ^ op (l,...,(a + 6r,...,fi) 

ph. pol. a 

+ Y SplitLt° P (a A °,^) A 1 ^(l,...,(a + br,...,' l 

ph. pol. (T 

We must again consider the singular channel specially. The cut in that channel is, 
Al-^(l,...,aA...,n)\ Sab cut = 

J dLIPS 4 - 2e (^,£ 2 ) A^ °v{-t u b + 1, . . . , a - 1, -^ 2 )4 ree (£ 2 , a, Mi) 

+ y dLIPS 4 " 26 ^!, ^ 2 ) ^TMi, b + 1, • • • , a - 1, -£ 2 )Al- loop (£ 2 , a, b, 4) 

+ y dups 4 - 2e (h,e 2 ,e 3 ) A^i-e^b + i,..., a -i,-e 2 , -e 3 )A^ ee (£ 3 ,e 2 , a , b, h) 

Once again, momentum conservation forces all legs crossing the cut to become collinear. For 
the two particle cuts, this is just the statement that — > 0, with the 4 massless. For the three- 
particle cuts, the three-particle invariant vanishes, tg 1 e 2 g 3 — > 0. Because the 4 are physical (in a 
D-dimensional sense), they all have positive energies; the vanishing of the three-particle invariant 
thus forces the separate two-particle invariants to vanish as well, se i e j — > 0. 

The behavior of tree amplitudes when three color- adjacent external particles become collinear 
is governed by analogs of the collinear splitting amplitude, 

*tree / \ a ,\ b \ c \ k a -k b , k a -k c ,k b -k c ^0 

A n+2 {...,a ,o , c ...) ► 

Y Split^V", &\ c^)A^{. . . , (a + b + of, . . .) , (5 - 4) 

ph. pol. a 



(5.3) 
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valid when n > 4. The squares of such double- collinear splitting amplitudes have been computed 
by Glover and Campbell [21] and by Catani and Grazzini [22]. 

As in the one-loop case, we cannot factorize the four- or five-point amplitudes on the right-hand 
side of the cut. We can, however, factorize the amplitudes on the left-hand side, obtaining 
£ 4^((a + &)*,& + l,...,a-l) 



ph. pol. 
0"Al > ^2 



x J dLIPS 4 - 2e (^,£ 2 ) Split^(-^-^)4 ree (£ 2 ,a,Mi) 
£ iCi((a + 6r,6 + l,...,a-l) 



ph. pol. 



x y dLIPS 4 - 2e (^,£ 2 ) SpUt^(-^-^)Af"(4,a,Mi) 
£ ^r i ((a + 6r,6 + l,...,a-l) 



(5.5) 



ph. pol. 



x J dLIPS 4 - 2e (^,£ 2 ) Split t _ r : e (-4 2 ,-^ 1 )^- loop (£ 2 ,a,6,£ 1 ) 
+ J] <- e i((« + &r^ + l,...,a-l) 



ph. pol. 
CT i A2 A3 



x y dLIPS 4 - 2e (£!,£ 2 ,4) Split^C-^ 2 , -4 3 , -^ 1 )A* ree (f 3 ,^2, a, • 

We recognize 

y dLIPS 4 - 2£ (^,^ 2 ) SpKt t ??(-e$ a ,-e$ 1 )AT e (t2,a,b,e 1 ) (5.6) 
as the cut of the one-loop splitting function, so that eqn. (5.3) can be written as 



J2 Split^V",^) A^((a + by,b+l,...,a-l) 

+ ^ Spht 2 _- loop (a A ",^) ^((0 + 6)^,6+1,. ...a-l) 



ph. pol. 
f i^l A 2 



(5.7) 



•Sab CUt 



where I have introduced a new quantity, the two-loop splitting amplitude Split 2 loop . It is defined 
by its cuts, 

Split 2 .-"** (« A -,^)|,., out 



dLn»S 4 - ae (£i^ 2 ) Split^ oop , )^ ree (£ 2 , a, 6, ^0 



ph. pol 



+ 



+ Split^-^ 2 , ~t X i )^4" lo ° P (^2, a, 6, £x)' 

^ / dups 4 - 2 ^e 1 ,e 2 ,e 3 ) s P \\t t z e (-£ X2 ,-£ X3 ,-£ x 1 1 )Af ee (e 3 ,e 2 ,a,b,e 1 ). 

(5.8) 
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ph. pol. 
ff ^li^2 A3 



The splitting amplitude has no cuts in other channels. Because the cut in this channel is taken to 
all orders in e, this equation suffices to reconstruct it completely. 




Figure 5. A schematic depiction of the collinear factorization of two- loop amplitudes. 

The first term in eqn. (5.7) repoduces the second term on the last line of eqn. (5.2); combining 
all channels, we finally obtain the collinear behavior of a leading-color two-loop amplitude, also 
shown schematically in fig. 5, 

A 2 n - loop (l,...,a K ,b Xb ,...,n) 

[spiiOa*- . ^) ^-T op u, •••,(« + &r, •••,") 

ph. pol. a 

+ Spht^ OOP (a A » , ^)C(1> • • • > (a + b)°, ■ ■ ■ , n) (5 " 9) 
+ Spht 2 _- loop (o A " , b^ ) AT-\ (1, ■ ■ ■ , (a + by , . . . , n) . 

6. One-Loop Corrections to the Double-Collinear Splitting Amplitude 

As three external legs become simultaneously collinear in a one-loop amplitude, we expect a 
factorization analogous to that at tree level, eqn. (5.4). We can derive it using an approach similar 
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to that in previous sections. Much of the discussion will carry over, but there are certain differences 
worthy of closer examination. 

The major difference is that the double-collinear splitting amplitude is no longer a function of 
a sole invariant; rather, it depends on the two neighboring two-particle invariants s a b and Sb c , as 
well as the overall three-particle invariant t a t, c . All of these invariants vanish in the limit. (I take 
the order of the external legs to be a, b, c.) As discussed by Campbell and Glover [21], the tree-level 
double-collinear splitting amplitude contains all terms that are singular in two of the invariants. 
(Terms singular in only one invariant would eventually give rise to a vanishing contribution to 
physical processes.) At one loop, we will acquire additional epsilonic powers of the invariants (and 
possibly more complicated analytic functions of their ratios) , but the basic requirement of having 
integral or half-integral power singularities in two of the invariants remains. 

As in the simple splitting amplitudes, these integral powers can arise either from amplitudes on 
either side of the cut, or else from the loop integral whose cut we are taking. In the former case, all 
three of the collinear legs must be on the same side of the cut. In the latter case, an inverse power 
in a given invariant also implies the integral function must contain a cut in that invariant. Let us 
focus on the inverse powers of t a b c - Two of the possible pole combinations — in (s a b,t a b c ) and in 
(sbc,t a bc) — manifestly require that the loop integral have a cut in t a bc- The third combination 
(s a b, Sbc) also requires such a cut. To see this, note that it must have cuts in both s a & and Sb c , and 
thus must have uncancelled propagators between legs a — 1 and a, and likewise c and c+ 1. These 
two uncancelled propagators force the integral to have a cut in the three-particle channel as well. 

The set of cuts we must consider is thus exactly analogous to that for the simple collinear 
splitting amplitude: all legs must be on one side of the cut. Contributions missed by ignoring cuts 
that separate the collinear legs will be be picked up in the primary singular channel, where the 
three legs will sit alone on one side of the cut. The remaining arguments about the order of limits 
go through as before: in channels other than the singular channel, we can take the singular limit 
before performing the cut, in the singular channel we must tread carefully. (These arguments also 
carry over to general multi-collinear splitting amplitudes.) 

In a non-singular channel, the limit of eqn. (3.1) again produces a tree-level splitting amplitude, 



> 

cut 



A^(l,...,a,b,c,...,n)\ td f 

^2 SpHt t Z e (a x ",b x >>,c x c) j 'dLIPS 4 - 2e (£i,-£ 2 ) 

ph. pol. <r • (QA) 

x AT- m+2 (h, d, ...,/, -£ 2 )AT(h, f+l,...,(a + b + cy,...,d-l, -h) 
= £ Split^,^) #(1 (fl + i + c)' n) 

ph. pol. a " CUt 
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In the three-particle singular channel, the cut is 
A^(l,...,a,b,c,...,n)\ tabc cut = 

E J dUPS 4 - 2 ^l 1 J 2 )ATJ > A-^ X \c+l,...,a-l,-£^)A^ 

ph. pol.Ai 2 

(6.2) 

The vanishing of the three-particle invariant forces the vanishing of the invariant of the legs crossing 
the cut, s^ 2 , and we can again use this to factorize the amplitude on the left-hand side, so the cut 
becomes 

E AT- e 2 ((a + b + cy,c+l,. ..,«-!) 



ph. pol. 



j dLIPS 4 - 2e (^,£ 2 ) Split^i-^*^-^ 1 ) Af 00 ^ 2 ,a x « ,b Xb ,c x %£^) (6.3) 



tabc CUt 



= £ <- 2 ((a + & + c)V+l,...,a-l) Split^ °>\^,c A 

ph. pol. a 

where the latter equation defines the one-loop double splitting amplitude. 

Although this splitting amplitude will also have cuts in the two-particle invariants within the 
collinear set, as discussed above it must have a cut in the three-particle channel, and hence this 
cut suffices to reconstruct it fully, 

split it OOI V a ,6 A v Ac ) = 

-A 2 _D— Ai^ 



. . , . ., . {£ + k a + k b + k c f 

ph. pol. Ai,A 2 

x Al ree {{-£ -a-b) X2 ,a x % b Xb , c Ac , £ Xl ) 

(6.4) 

The restriction to physical polarizations will once again give rise to transverse projection operators 
inside the loop. 

Combining all cuts, we find a double-collinear factorization formula very similar to the simple 
collinear one, 

(X. II foil c 



^- loop (l,...,a A »,6 Ab ,c Ac ...,n) 



£ (spht^V , &\ ^-° 2 op (l, • • • , (a + 6 + c) CT , . 



,n 



(6.5) 



ph. pol. a 

+ Split ^°° P (a A " , ^,c A ^I r ! 2 (l, . . . , (a + b + c) CT , . . . , n 
7. Collinear Factorization to All Orders 

In this section, I give a proof of collinear factorization to all orders, generalizing further the 
results of the previous sections. As in the two-loop case, I will restrict attention to leading-color 
amplitudes, that is to coefficients of N l times the tree-level color structure in an Moop amplitude. 
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The general cut formula here is 
A l - lo °v(l,...,a u a 2 ,...,a c ,...,n)\ td f 

l + l n 

Y / dLIPS 4 " 26 ^!,...,^) 



3=2' 

X Y A n-ZM -,f,tl,-,ti) ^+ 1 r'" fcH ° OP (/ + l, '/ 1- 

k=0 

(7.1) 

so long as the dj all end up on one side of the cut. (As in earlier sections, we can ignore contributions 
where different cij end up on opposite sides of the cut.) The outer sum is on j-particle cuts, the 
inner one on the different j-particle cuts; m is as in previous sections, the number of external 
momenta on the far side of the cut; and '0-loop' means 'tree'. The sum over intermediate particle 
types and physical polarizations is, as always, implicit. 

The proof will proceed by induction, in fact a double induction, over both the number of loops 
and the number of color- adjacent legs becoming collinear simultaneously. Denote by S c the set 
of invariants built out of consecutive momenta in {tti, . . . ,a c }, and indicate the multiple collinear 
limit by S c — > (that is, all of the invariants are supposed to vanish). 

Now assume that for r < I and c < n — 2, that we have the following factorization, 



A n p (l, • • • , a>i 1 , & 2 2 , • • • , a c ° , . . . , n) 



> 



Y £ S P hr_t op (a^ , a 2 A2 , . . . , «£" ) 4-" c t loop (l, . . . , (ax + . . . + a c ) CT , . . . , n) , (7 ' 2) 

ph. pol. a v=0 

The starting point for the induction — r = 1 and c = 2 — is the one-loop factorization proven 
in section 4. In addition, I will make use of the tree-level multi-collinear analog to eqn. (5.4), 
which can again be derived either from the Berends-Giele recurrence relation or from the string 
representation. 

In all cuts except that in the primary singular channel = (k ai + . . . + k ac ) 2 ), we then find 
using the above assumption, that 



/W-1°°PM n x i r M 

j± n ^i, . . . , a 1 , u 2 , • • • , 

l+i , l+i-j 



Ac „\\ S c ^0^ 

t d ...f cut 



a c ,...,n)\ 



Y / dLIPS 4 - 2 ^,...,^) Yl f,*u-,tj) 

3=2 fc=0 (7 _ 3) 

x Y E spht!:^^ 1 ,^,---,^) 

ph. pol. a v=0 

x AilYJ c - k - v) - loop (f+l, . . . , ( 0l + . . . + a c r, ....,/ 1. t 3 -h) . 
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Interchanging the summations, 

l+i l+l-j l+l-j-k l+l l+l-j l+l-j-v l-l l+l-v l+l-j-v 

EE E=EE E=EE E . <w 

3=2 k=0 v=0 j=2 v=0 k=0 v=0 j=2 k=0 

we can transform eqn. (7.3) to obtain 
A l - lo °v(l,...,a^,a¥,...,a^,...,n) 



td-.-f cut 

l-l l-v+1 



E E Split_t° P (a* 1 , a 2 A V . . , a A < ) £ / ^LIPS 4 "^ (*,..., ^ 



ph. pol. £7 f=0 j=2 

x ^+ 1 .^; fc r^- loop (/+l, . . . , (ai + . . . + a c ) CT , ....,/ 1. 

= E E S P lit - OP ( a i 1 ' £t 2 2 ,---,«c Ac ) ^::]_f op (l,...,(a 1 + ... + a c r,...,n) 

— — td---f cut 

ph. pol. cr u=0 

(7.5) 

which is exactly eqn. (7.2) in this channel, but now for r = I. It again holds for all c < n — 2. 
In the primary singular channel, we obtain in the limit, 
l+i f . i+i-j _ k 

E / dLiPs 4 - 2e (^,...,^) e E E s P nt - oop ^'---'^) 

3=2 J k=0 ph. pol. av=0 

x A^- C ti 00p (a c +1, . . . , ai -l, + . . . + a c ) CT ) 

X ^c+j J ^ P ( a l 1 i a 2 2 ' • • • i a c% ~~ ^ji • • • i ~ ^l) 

= E E E4 fe -cS l0 °>c+l,...,a 1 -l,(a-l + ... + a c n 

fc=0 ph. pol. a v=0 
l+l-k r 

x E / ^LIPS 4_2e {h,...Jj) Split^°° P {l U ...,lj) 

j=2 J 

X -^c+j J P ( a l 1 > a 2 2 J • • • i a c C > ~~ ^j) • • • ) — ^l) 

which gives us an explicit formula allowing the reconstruction of the /-loop multiparticle splitting 
amplitude in eqn. (7.2) as the coefficient of A^ c+1 , 

spiit^r^,^...,^)! 

\Kj. cut 
l-l l+l-k . 

= E E E y dLips 4 - 2£ ^i, . . . ,t s ) spht fc _r p (^r CTi , . . .,ej aj ) ( 7 - 7 ) 

fc = .7 = 2 ph. pol. &i 

v 4('+l-j-fc)-loop/ Ax A 2 A c 

A ^l c+ j V"l ) "2 5 • • • > U C 5 * j 1 ■ ■ ■ 1 t l ) 

Combining the two kinds of channels, we obtain eqn. (7.2) uniformly for r = I, thereby proving 
the desired result. 



(7.6) 
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8. Conclusions 



The square-root factorization of amplitudes in the collinear limit is a striking feature of gauge 
theories, yet one that hides subtleties. The subtleties, and the complexity of a conventional di- 
agrammatic approach to the problem, are associated with the presence of infrared singularities. 
The naive intuition of computing only those diagrams associated, for example, with fig. 2, is not 
quite right. This is demonstrated graphically by the presence of dilogarithms in the quark- gluon 
splitting amplitudes [7,8], which cannot arise from massless three-point integrals at one loop. The 
unitarity-based proof given in this paper circumvents these complexities, and gives as well explicit 
formulae for computing the splitting amplitudes. 

I thank P. Uwer, Z. Bern, and L. Dixon for helpful comments. 
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